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1 Introduction 

The "final fate" of gravitational collapse of spherically symmetric dust clouds is well under- 
stood. The central singularity can be naked or covered depending on the choice of initial 
data; there is a "critical branch" of solutions where a transition from naked singularities to 
black holes occurs (see [||, § and references therein). 

Recently, a considerable effort has been paid in order to extend the above mentioned 
results to collapse models with general matter fields. In particular, one would like to under- 
stand the role played by stresses in the formation and the nature of the singularities, since 
both numerical 0, 11 1^ and analytical P, ^ studies tend to confirm violations of the cosmic 
censorship conjecture of the kind observed in dust spacetimes. 

Gravitational collapse with stresses is quite difficult to handle analytically even in spher- 
ical symmetry, one of the reasons being the lack of physically valid exact solutions. In a 
recent series of papers P, Q we have obtained the general exact solution for the case of 
anisotropic materials sustained only by tangential stresses (see |l^ for a review on the role 



of anisotropic materials in relativistic astrophysics); the nature of the singularities forming 
in such solutions is still largely unknown (exceptions are the self-similar case and models 
with special equation of state [|11|, 0). 

Among the solutions with vanishing radial stresses, there is a system of counter-rotating 
particles - the "Einstein cluster" - first introduced by Einstein |]T3 (see also W^) in the static 



case and then generalized to the non-static case by Datta [T^ and Bondi |jT6| (see also JT? 



and [|TB]). The motion of the particles in the cluster is sustained by angular momentum whose 
average effect is to introduce a non-vanishing tangential pressure in the energy-momentum 
tensor. This model is particularly interesting from the physical point of view since it can give 
some hint on the effects of rotation on collapse, without raising difficulties usually associated 
with the axisymmetric spacetimes. The analysis of the singularities of the Einstein cluster 
model has been recently pioneered by Harada, Iguchi and Nakao ( |]19[, to be referred after- 
wards as HIN). Using the "mass-area" coordinate technique, developed in 0, these authors 
discovered a particular solution within the Einstein cluster class for which the quadratures 
can be explicitly carried out. This allowed them to discuss the nature of the singularity 
for this special solution in full details. They have also discovered several features of the 
singularities forming in the general model, for which the metric is not expressible in terms 
of elementary functions. Their results comprise several new examples of naked singularities. 
Our results here extend the HIN results to non-analytic cases. We show that the initial 
data can be parameterized in terms of the first non-vanishing derivatives near the centre, 
as in the dust case, and obtain the corresponding spectrum of endstates of the collapse 
for marginally bound initial configurations. Several new interesting features of the Einstein 
cluster model arise. 



2 The Einstein cluster model 

The general solution of the Einstein field equations for spherically symmetric collapse with 
tangential stresses can be reduced to quadratures using a technique first introduced by 



Ori pO| to study the spherically symmetric charged dust. We will very briefly recall the 



solution here, referring the reader to paper ^ for details. 

Systems with vanishing radial stresses are characterized by conservation of the Misner- 
Sharp mass (m). Therefore, one can use it as a comoving label and introduce a system of 
coordinates ( "mass-area coordinates" ) in which the line element turns out to be 
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Here K, u and h are functions of R and m. The function u is the modulus of the velocity of 
the collapsing shells and satisfies 



while 



where 



, 2m ^2 



K = g{m)+ f G{m,R)dR 



G{m,R) :-- 



RE 



1 + 



R 



E^ 



2m 
~R^ 



E^ 

/l2 



-3/2 



(2) 



(3) 



(4) 



The functions g{m,) and E{m) are arbitrary. The quantity /i(m, R) is a measure of the 
internal energy stored in the material, and therefore plays the role of equation of state. It 
can be shown that the quantities characterizing a material with vanishing radial stresses, 
namely the energy density e and the tangential stress 11, can be written in terms of h as 

follows m- 



h 



AttuKER^ 



j.__Rdh 
~~2hdR^ 



(5) 



These formulae show that the tangential stress vanishes whenever the function h does not 
depend on R. Therefore in this case the material is not sustained by any stress and as a 
consequence the line element (|lD reduces to the dust one (Lemaitre-Tolman-Bondi space- 
time). Since, in all formulae, E and h appear only in ratio, the value of h on an arbitrary 
curve R = Ro(m) can be rescaled to unity, so that, in particular, the dust solutions can be 
characterized by /i = 1. The function Ro{m,) plays the same role as played by the initial 
mass distribution in the "standard" comoving coordinates (the inverse transformation from 
mass label to lagrangian label being r = RQ{m)) while the function E{m) plays the role of 
"energy function" familiar from dust models (see for details). 

The above recalled structure shows that a solution with tangential stresses is identified 
(obviously modulo gauge transformations) by a triplet of functions {g, E, h}. This "param- 
eterization" contains the dust spacetimes as the "subspace" {g,E,l}. In this way we can 
construct in a mathematically precise way, a comparison between dust and tangential stress 
evolutions by comparing the end state of the dust collapse {{g,E,l} with chosen g and E) 
with the endstates of the tangential stress solutions {g, E, h} with the same g and E and 
different choices of the equation of state h (see |^, ^ for details). 



The physical singularities of the spacetime described by the metric (|1]) correspond to 
infinite energy density and are given by-R = Oorii' = 0. At-R = the shells of matter 
"crush to zero size" (shell-focusing singularities) while vanishing of K implies intersection 
of different shells of matter (shell-crossing singularities). We are interested here only in the 
central shell-focusing singularity. It is important to distinguish the central singularity (i.e. 
that shell focusing singularity occurring at i? = m = 0) from the non-central singularities 
(singularities occurring at i? = 0, m 7^ 0) which are always covered in spacetimes with 
vanishing radial stresses, as a simple consequence of mass conservation. 

At i? = m = the apparent horizon and the singularity form simultaneously. Therefore, 
to analyse the causal nature of the singularity, one has to study outgoing radial null geodesies 



and check if there are some which meet the singularity in the past. This can be done p2 
by defining a quantity x = R/2rrf (where 1/3 < a < 1) and searching for a real positive 
solution xq of the following "root equation" 
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This equation depends both on the choice of the equation of state h and on the initial data. 
Setting h = 1 one can deduce from it the final fate of the dust gravitational collapse P, p3|] , 
so that it can be used to study the way in which the end state of collapse for fixed initial 
data is altered by a non-trivial h. 

To analyse the roots equation, one faces with the fact that the integral appearing in 
formulae (P) introduces a sort of non-locality which is difficult to handle. As a consequence, 
it is very difficult to work in full generality and only a few cases have been analysed so far. 
Among different cases with tangential stresses, a special role is played by a system of counter- 
rotating particles (The Einstein Cluster). This system has several attractive features from 



the physical point of view since it mimics in a simple way the effect of rotation (see e.g. ||24|| ). 
So far, the only known results on the nature of the singularities forming in this model have 
been obtained in the HIN paper (such results will be recalled in the next Section). 

It has been shown in paper P] that the Einstein cluster model can be obtained from the 
general exact solution with tangential stresses by choosing a specific form of h: 
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In the above formula, the functions Hq and L are arbitrary. The function L is the specific 
angular momentum of the particles while only one among the functions E and Hq is inde- 
pendent. Therefore the function h^ can be used to normalize the value of h to unity on the 
initial data surface R = Ro{m). This choice does not affect the generality of the equation of 
state (0) within the Einstein cluster class until the function E is chosen. 

The final fate of the collapse turns out to depend on the behaviour of the physical 
quantities near origin, on the initial data surface. Three such quantities are independent, 
and we shall use from now on F, L^ and / = E"^ — 1. We assume these functions to be 



Taylor-expandable with respect to r (the "standard" radial coordinate r = Rq) near r = 0. 
It is now easy to check (see |§]) that regularity at the initial surface requires the following 
kind of expansions: 

F(r) = For3 + F„r"+3 + ... (8) 

/(r) = // + ... (9) 

L\r) = Lur^ + ... (10) 

where I > 2, k > A and dots denote higher orders terms. The n-terni in F{r) has been put 
in evidence because it characterizes the first non-vanishing derivative of the initial density 
profile; a positive Fq and a negative F„ are required for physical reasonability. 

In the above formulae we allow for terms of any parity. It is very important to specify 
carefully the meaning of this choice, since a considerable debate on this point has already 
taken place in the case of dust (i.e. L^ = here). 



In the dust case, pioneering numerical results by Eardley and Smarr |25| were first inves- 



tigated analytically by Christodoulou |2^ and Newman [^. Both these authors required 
the density and metric functions to be smooth {C°°) with respect to a system of cartesian 
coordinates near the centre (a function of this kind is sometimes called extended analytic). 
Obviously, this implies that a power series expansion near r = can contain only even powers 
(for instance, if a linear term is absent but a cubic term present in the mass F{r), the initial 
density e{r,0) = F'/ATrr"^ is only C^ near the centre w.r. to a local cartesian system). From 
the physical point of view, however, this restriction is unnecessary and actually rules out 
many interesting situations [^ . Indeed, the results by Christodoulou and Newman were 
extended to the case of an arbitrary parity by Singh and Joshi ||2^ showing analytically, for 
instance, the existence of a transition parameter at ra = 3 for marginally bound data (this 
structure will be very briefly recalled in the next section). 

Returning to the case of Einstein cluster, the extended analytic case here corresponds to 
even values of k in the expansion of L^ and was first considered in the HIN paper, while we 
allow here for odd values of k as well as for even ones. 

The expansions can now be translated in mass-area coordinates. In particular, we put 

L\m) = P^m'^/' + ... , (11) 

where Pk is a positive constant. 



3 The endstate of collapse 

The dust models admit no globally regular solutions. This is not true for general matter 
fields even within spherical symmetry. For example, in homogeneous perfect fluid models, 
we can have globally regular solutions which satisfy suitable energy conditions and have 
regular initial data (see e.g. |^); in other words, the endstate of a non-static model can 
be an eternally oscillating solution or a bounce-back scenario. It is, therefore, important to 
understand the qualitative behaviour of the motion as a preliminary step before investigating 
the causal structure. 



The equation of motion (^ governing the collapsing shells can be conveniently written 
as 

2 Z{R,R^) 
^ R{R^ + L^)' ^ ' 

where the "effective potential" Z is defined by 

Z{R, Ro) = fR^ + 2mi?2 + L2(2m - R) . (13) 

This function should be interpreted as an analogue of the Newtonian effective potential gov- 
erning the motion of the fixed shell Ro, so that the allowed regions of the motion correspond 
to Z > 0. It is immediately seen that the sign of / = E'^h'^ — 1 governs the behaviour at 
large values of R, since unboundedly large values of R are permitted only for / > 0. The 
region near i? = is, instead, always allowed for any Rq ^ since Z{0, Rq) = 2mL'^ > 0. 
To investigate on the formation of the central singularity (i.e. the singularity forming at 
R = Rq = 0), let R = ((t, RqjRq. To the lowest order in Rq one has 

Z ^ i-MRo'-' + 2Foe + foC + 2P,FoRo'-' + PkC'Ro'-' + /3fc/oC'i?o'-')i?^ , (14) 

(recall that F = F^R^, f = /o-Rg ^^^ -^^ — PkRo to the lowest order). Near C = 0, the above 
equation gives Z ^ —P^CRq < for fc = 4, and therefore the singularity does not form. 
Since it is the positivity of /34 which does not allow singularity formation, tangential stresses 
regularize the corresponding (/54 = 0) dust solutions. For k > 4 one has Z ^ 2Fq('^RI > 
and the singularity always forms (the case k = 4 was firstly studied by Evans ||T^ , while the 
HIN paper contains the analytic case of even k). 

For A; > 4 we have to analyse further the dynamics of the non-central shells, which 
depends on the existence of roots of the cubic equation Z = 0. Due to the above remarks on 
the behaviour of Z for big and small values of R, we can a priori expect the following kinds 
of situations. If Ri, R2 and i?3 are the roots of the cubic and Ri < R2 < R3, one has the 
following sets of allowed regions 

0<R<Ri, R2<R<R3, / < 0, 

0<R<Ri, R2<R, / > 0, (15) 

where, of course, the forbidden regions between R2 and R^ can disappear if either one or 
both of these two roots vanish or are not positive (thereby giving < R < Ri for / < and 
i? > for / > 0, respectively). 

The fine details are quite involved and were given for the first time by Bondi |T^. The 
point of view here, however, is slightly different since we are interested only in the dynamics 
of the shells near the central one to check whether the central singularity becomes trapped 
or not. A simple way to do this is to evaluate the effective potential on curves of the type 
R = XF with A > 2. For the region near centre we obtain 

Z{XF, Ro) ^ (2 - X)f3kFoR'^+^ + 2X'F^RI (16) 

For k < 5 the first term is the leading one and it is negative. For A; = 6 it is possible to find 
values XF such that Z{XF, Rq) is negative if the quadratic equation (2 — X)PqFo + 2FqX^ = 



has a positive root. This in turn imphes that the quantity D := yf^/ Fq must be greater 
than four (HIN). In both these cases (fc = 5 and k = Q for D > 4) the central singularity 
is naked. Indeed, \F is not an allowed value since Z is negative at it. But F goes to zero 
as i?0) ^^d therefore the shells near the centre must have initial data R = Rq greater than 
2F. For sufficiently small Rq, this will hold true also for XF, so that a region of the kind 
R > R2 > 2F containing R = Rq must exist. But the corresponding dynamics can never 
cross the forbidden region to reach values smaller than 2F. Therefore, the region near the 
centre is untrapped eternally in these cases: the apparent horizon does not form and therefore 
the singularity is naked. This phenomenon has been firstly discovered in the HIN paper for 
the case k = 6,D > 4. Interestingly enough, this kind of eternal singularity is completely 
different from the naked singularities arising in dust spacetimes, which are usually visible 
only for a finite amount of time before apparent horizon formation. A peculiar characteristic 
of the singularities here is that they are not accompanied by non-central singularities, since 
all the shells near the central one remain regular eternally (the region near i? = is forbidden 
for them). 

For k = 6,D < A as well as for any k > 6 the above argument cannot be applied. The 
special case k = 6, D = A, f = can, however, be treated analytically since the integral 
appearing in formula (|^) can be expressed in terms of arbitrary functions. The singularity 
in this case turns out to be always naked; again independently from the choice of the initial 
matter distribution (HIN). 

To understand the nature of the singularities in the remaining possible cases, we have 
to resort to the root equation. The integral which appears in formula (^ cannot be cal- 
culated explicitly, however, the existence of a naked singularity can be probed as follows 
(for simplicity we shall consider the case E{m) = 1). We expand the integrand near centre 
in the powers of m and express the integral as the sum over successive integrated terms of 
the series. The proof that this is possible makes use of Lebesgue's dominated convergence 
theorem, as shown in details in the Appendix. 

For k = 6 equation ([A6|) gives 
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where dots stand for terms going as 0{m~^^\ In this equation it is possible to fix uniquely 
the value of a for each choice of n, in such a way that the resulting algebraic equation in 
Xo has a positive root (see Table 1). Therefore the singularity is always naked for k = 6, 
independently from the choice of the initial density distribution. 

As recalled in the previous section, the value of the root gives the tangent of the escaping 
geodesies near the singularity. It is somewhat unclear if a physical content can be given 
to the actual ("numerical") value of this tangent, although it is very likely that such an 
interpretation could be found in the near future. In any case, we can get some insight into 
the effect of the state equation on naked singularity formation looking qualitatively to the 
role of the parameter Pk in the roots. With this aim, the values of the roots for the various 
possible initial distributions of mass in presence of a non-zero Pq (counter-rotating particles) 
can be compared with the corresponding values for dust (/Jg = 0) (see Table 1). 
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Table 1: The endstate of counter-rotating particles for A; = 6 vs dust 

At n = 1 the inhomogeneity dominates over the effect of counter-rotation: both the 
endstates are naked singularities and the roots coincide. At n = 2 a contribution coming 
from the counter-rotation appears and both the endstates are still naked singularities. At 
n = 3 the root equation for dust spacetimes becomes a quartic which has positive roots 
only if the a-dimensional quantity ^ := F3/{2^/^F^^^) is less than ^^ = -(26 + 15y3)/2 Q. 
Therefore, C,c can be called critical parameter: at ^ = ^c a transition occurs and the endstate 
of collapse turns from a naked singularity to a black hole. In the case of the Einstein cluster 
instead the root becomes independent from n at n = 3 and is always present for /^g > 
making the endstate always a naked singularity. For n > 3 dust always forms a black hole 
while here the endstate is always a naked singularity. 

The existence of such "tenacious" naked singularities which are not sensible to changes 
in the initial mass distribution, first discovered in the HIN paper, raises the question if these 
singularities are already present in dust models or are peculiar of spacetimes with stresses. 
Actually it turns out that such special dust spacetimes exist, as we are going to show. 

Observe that the energy function of the HIN model is of the form /(r) = /4r^ -|- ... (where 
fi = — 2^Fq ), i.e., the expansion of this function starts from the quartic terms. Now one can 
ask for the final fate of those dust spacetimes which have initial velocity profiles of this kind. 
To the best of our knowledge, the answer to this question is not present in the literature since 
all previous studies on dust spacetimes have assumed a non-vanishing second order term. 
The causal structure can, however, be analysed easily and it turns out that such solutions 



actually have the same causal structure of HIN: the singularity is always visible |^| 
It remains to consider the cases k > 7. For k = 7 formula (|A6|) gives 
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where dots stand for terms going to zero as 0{m 



k+n-7 , 



This equation leads to a structure 



of the endstates which is similar to that occurring in dust (second column of table one) 
for n = 1 and 2. At the transition {n = 3), the region of naked singularities formation 
is characterized by the inequahty ^ < ^c where ^ = ^ — 3/37/(5-\/2Fq ) and ^ is the dust 
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parameter recalled above. Since [3i is positive, a sector of the black hole region in dust is 
uncovered by the effect of counter-rotation; for n > 3 it remains uncovered for strong enough 
rotation parameter {(3-j/Fq > — 5v2^c/3). 

Finally, for k > 7 the causal structure of the singularity is identical to that of marginally 
bound Lemaitre-Tolman-Bondi spacetimes, since all the terms coming from rotation vanish 
near the singularity (a similar situation has been recently found to occur in a completely 
different model with tangential stresses [0]). 



4 Discussion 

In recent years, a big effort has been undertaken to understand the final fate of the grav- 
itational collapse of a dust cloud. Indeed, one can safely assert that this end state is now 
completely known in dependence of the choices of the initial data. 

The dust model is the most simple model which can be conceived in gravitational collapse 
and in the absence of general proofs (and indeed even of mathematical formulations of) of 
cosmic censorship theorems, it was meant to be used as a tool to get insights in to more 
general collapse situations. We are, therefore, in some sense ready to approach the problem 
with general stresses and indeed some analytical results are beginning to be known. In this 
respect our full understanding of dust spacetimes should hopefully be used as a starting 
point. 

As a simple generalization of the dust models, we proposed to analyse spacetimes with 
tangential stresses and constructed a tool for doing that using mass-area coordinates. The 
first results on the nature of the singularities for such spacetimes were obtained in the HIN 
paper for the case of the Einstein cluster model. In the present paper we have extended the 
HIN results to the non-extended analytic case. 

The picture which arises is very intriguing. The key role is played by the parameter k 
measuring the "strength" of angular momentum near the centre (L^ ^ A-Rq ^^ -^o tends to 
zero). For k = 4 the solution is regular and tangential stresses avoid singularity formation. 
For fc = 5 we have central singularity formation, as in dust, but the neighbouring shells are 
regular; all such shells are untrapped forever allowing the central singularity to be eternally 
visible. This phenomenon is clearly connected to singularity theorems since all the solutions 
considered here satisfy the energy conditions and therefore, in the presence of trapping the 
shells, would have to become singular. For /c = 6 a sort of transition takes place, since for 
D = \/j3e/Fo < 4 also non-central shells can become singular, the exact HIN solution playing 
the role of separation point. For k = 7 the remaining effect of counter- rotation becomes very 
weak and it is effective only at the n = 3 transition from black holes to naked singularities, 
shifting the transition parameter C,c and uncovering part of dust black hole region. The naked 
singularities continue to exist for strong enough rotation parameter /3j. For k > 7 the final 
fate is the same as in dust models. 

From the physical point of view, this picture is due to the non-convexity of the state 
function which acts as a "source of acceleration" near the centre avoiding the formation 
of apparent horizon and therefore enforcing nakedness. The structure of the spectrum of 
endstates thus depends in a peculiar way on the fact that the Einstein cluster model, when 
interpreted as a material characterized by a macroscopic equation of state, does not admit 

9 



a local minimum of the energy density. This is readily seen since the derivative of the state 
function h with respect to R is strictly negative and proportional to —L^ near the centre, 
while equations of state which have to be expected for, say, nuclear matter at high densities 
should admit a status of (local) minimization of internal energy. However, the HIN results, 
as well as our results here, do raise the interesting question if the effects of rotation might 
become dominant with respect to those of a non-trivial equation of state at high densities; 
when the hypothesis of average geodesic motion (which is the basis of the relative simplicity 
of the counter- rotating particles model) is lost. To approach this question one should turn 
to the (so far unsolved) difficulties of the non-spherical collapse scenario. 

Remaining in the safety zone of spherical symmetry, our results here support the view 
that formation of naked singularities has to be expected in generic situations of spherical 
collapse [0, 0. 
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Consider the function 
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where G{m,R) is defined by equation (^). Due to the Mean value theorem there exists 
x(^w) G {RQ,2m°'x) such that 
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Since Rq ~ m^^^ as m goes to zero, 1/3 < a < 1 and x is finite and positive, both m ^/^x(m) 
and m°'/x{'m) are finite and positive as m goes to zero. Therefore we can evaluate the right 



hand side of equation (|A2|) using (|J) as follows 



-Ro 



2^2(1 + /3fcm('=-2)/3Fo/3^i/2 V m^^^ 



l + Z^.^e^-e-Vaf^ 
X{m) 



+ ^^(^-4)73 fxM^ p2/3 _ Pk^(k-3-3a)/3 



X 



1 + /3km 




m^ 



m 



Xim) 



+ 



-3/2 



(fc-6Q)/3 / ^ 



+ ^(A;-2)m 
6 



(fc-4)/3 f xim) \ 2/3 

ml/3 i ^0 



A 



fc ^(2A;-5-3a)/3 



m 



xim) 



p2/3 _ k(3^ (k~3-3a)/3 

° 6 



m 



X{m) 



+ 



(A3) 



10 



where dots stand for terms of higher order in positive powers of m^/^. For A; > 6 only the 
factor m~^ is divergent, all the other terms being finite due to a < 1. Therefore, the limit 



/2(wi; a) := lim 



m^rO 



2m°'x 
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is convergent. Hence, using Lebesgue's dominated convergence theorem, we can expand the 
integrand near the centre {m = 0) in the leading powers of m and integrate successive terms. 
This gives 
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Using this equation in (^), we can write the root equation in the form (valid for k>Q) 
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